Boundary spin Hall effect in a two-dimensional semiconductor system with Rashba 

spin-orbit coupling 
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We derive boundary conditions for the coupled spin-charge diffusion equations at a transmitting 
interface between two-dimensional electron systems with different strengths of the Rashba spin- 
orbit (SO) coupling q, and an electric field parallel to the interface. We consider the limit where 
the spin-diffusion length l 3 is long compared to the electron mean free path I, and assume that a 
changes discontinuously on the scale of l s . We find that the spin density is also discontinuous on 
the scale of l s - In the case where the electron mobility is constant across the interface, this leads 
to the complete suppression of the expected spin injection from a region with a / into a non-SO 
region with a = 0. 

PACS numbers: 72.25.-b,72.25.Mk,73.23.-b,73.50.Bk 



I. INTRODUCTION 

Developing robust mechanisms for spin generation and 
detection is a central objective in the rapidly growing 
field of spintronics i There are many exciting venues 
for basic physics and nanotechnology to exploit in- 
jected spins in various systems, in both the semiclas- 
sical transport regime and the quantum regime where 
spins can encode entangled quantum bits. Since the 
early days of spintronics, 2 however, exciting theoreti- 
cal predictions often seemed to be plagued by experi- 
mental ambiguities. A notable exception is the field of 
metallic magnetoelectronicsj 3 - where many of key exper- 
imental findings have been well understood and the the- 
oretical predictions in turn were instrumental in push- 
ing the experimental frontier. Semiconductor spintronics 
which possesses much richer phenomenology and appears 
more desirable for technological applications, however, 
has caused more problems for theorists and experimen- 
talists alike. 

In systems with intrinsic spin-orbit (SO) coupling, the 
spin degrees of freedom are intricately entangled with 
the orbital motion, and even the concept of spin cur- 
rent (defined as a symmetrized product of spin and ve- 
locity operators) has resisted the thorough qualitative 
understanding^ There is an opinion that for many ex- 
perimental implications, it is more relevant to calculate 
spin densities rather than spin currents. This is a princi- 
pal motivation of this paper. The spin-density generation 
and dynamics in diffuse bulk semiconductors are conve- 



niently described by the semiclassical diffusion equation 
which is usually derived in the limit of weak SO coupling 
and dilute disorder .Si&iL&ui A, r _1 <C Ep, where A is the 
characteristic SO splitting, r _1 is the impurity scatter- 
ing rate, and Ep is the Fermi energy (setting h — 1 here 
and henceforth) . In order to calculate the nonequilibrium 
spin accumulation generated by the spin Hall effect at a 
Hall-bar edge or at a boundary between two different con- 
ductors, the diffusion equation has to be supplemented by 
the appropriate boundary conditions (BC) for spin and 
charge densities. The latter have been the source of a 
vigorous discussion in recent literature^^ ' 10 i 11 : 12 In ad- 
dition, the boundary problem with SO interactions was 
previously discussed in the context of spin-polarized bal- 
listic beam reflection and refraction 13 and in the regime 
of Friedel-like spin-density oscillations on the scale of the 
Fermi wavelength at sharp boundaries^ 

In this paper, we develop a general approach for sys- 
tematically deriving semiclassical BC for diffuse spin 
transport with weak intrinsic SO interaction. We assume 
that EpT ^> I, but we work in the "SO dirty limit," 
where Ar <C 1 and the D'yakonov-Perel spin-diffusion 
length l s is much larger than the electron mean free path 
I. (This limit is, in fact, necessary for the diffusion equa- 
tion itself to be meaningful near a boundary.) We use 
a Keldysh kinetic equation approach, 15 similar in spirit 
to the one used in Ref. [H for deriving the spin-diffusion 
equation. Our method reveals the physical meaning of 
the BC, which possesses an interesting correspondence 
with the classical Hall physics. We believe that our find- 
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FIG. 1: Simple example: A two-dimensional electron gas in 
the xy plane with a uniform electric field E in the x direction, 
and an inhomogeneous Rashba coupling a(y), which varies 
between a = ati, for y < 0, and a — ata, for y > d. If 
the in-plane spin polarization s y , which occurs in the region 
y > d spills over into the normal region y < 0, there will be a 
nonzero value of s z near the boundary, which can be measured 
optically. We find, however, that if the electron mobility is 
independent of y, there is no spillover of s y , and s z = 0. 



ings resolve some discrepancies existing in the literature. 



II. KINETIC EQUATION APPROACH 

Of central importance as a convenient case study as 
well as an experimentally relevant model is the noninter- 
acting Rashba Hamiltonian which describes the simplest 
two-dimensional system with intrinsic SO coupling: 12 
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H so = ~^{a{r),p x }a y + i {a'(r),p v } a x , (1) 

where p = p — eA, p = — iV is the canonical momen- 
tum, e the electron charge, A a vector potential, to the 
effective electron mass, a a vector of Pauli matrices, U 
the potential due to disorder and external gates, and a 
is the Rashba SO interaction parameter, which here de- 
pends on the position. The second coupling constant a' 



is equal to a for pure Rashba coupling, but the situation 
a' 7^ a applies, e.g., to structures on a [001] surface of 
GaAs, when linear Dresselhaus coupling is present, if the 
x axis is chosen along the (110) crystal direction. Note 
that we are disregarding the SO coupling due to the lat- 
eral potential U, and the electrons are moving in the xy 
plane. 



Although our approach is general and can be applied 
to other SO Hamiltonians and different heterostructures, 
we will focus on the BC in the Hall configuration of Fig.[T] 
Specifically, we assume a uniform electric field E = — <9t A 
parallel to the x axis, with a reflecting wall at some nega- 
tive y blocking any transverse charge currents. The two- 
dimensional electron gas is translationally invariant along 
the x axis, apart from the disorder potential, but the sys- 
tem parameters depend explicitly on y. We assume that 
the Rashba parameter a has different constant values ot\ 
and «o for y < and y > d, with a "transition region" 

< y < d. For reasons that will become clear later, we 
choose d to be long on the scale of the mean-free path 

1 = vpT [vp = ^/2Ep/m is the Fermi velocity), but short 
on the scale of the D'yakonov-Perel spin-diffusion length 
l s = 1/ A.qt (with Ao = aoPF and pp — mvp), so that 
the spins do not relax while traversing the transition re- 
gion. Note that by definition, l s = ^Dt s , where D is the 
diffusion constant and r s a characteristic spin-relaxation 
time. A geometry of this type can be produced by a 
charged gate covering one half of the 2D system, say the 
region y < 0, in which case the transition length d will 
be determined by the setback distance of the gate. In 
such a system, not only a, but also the potential U, the 
electron density n, and the electron mobility er/m will 
generally depend on y. However, for pedagogical pur- 
poses, we shall first consider a simplified model, where 
we assume that only a depends on y. Furthermore, we 
assume pure Rashba coupling, and we take a\ = 0. 



The kinetic equation with y-dependent a can be de- 
rived using the approach already discussed in Ref. with 
the result 
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containing an additional term proportional to d y a, which 
is at the heart of our discussion. Here, fj = z x <r, with 
z the unit vector along the z axis, and 



P{£) = to J (2^ 9(P ' 



(3) 



The second term on the left-hand side of Eq. Q is the 
drift with spin-dependent velocity, the third term is the 
spin precession in the Rashba field, and the fourth term 
is the spin-dependent acceleration in the y-dependent 
Rashba field. The right-hand side is the collision inte- 
gral in the self-consistent Born approximation due to the 
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weak and isotropic scalar disorder. A — i(G R —G A ) is the 
spectral function in terms of the retarded and advanced 
Green's functions G R ' A , which is given by 



A(p,e) 



no 



(4) 



In equilibrium 
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g (r,i;p,e) = - 



tanh(/3e/2) 



l(r,t;p,e) . 



(8) 



Here, 



A s (p,e) = 



1/t 



(l/2r)2 ' 



(5) 



Ve = V ' p/p = sm S <J X ~ cos <Tj,, where is the angle of 
p with respect to the x axis, and 



e ps = o E F + s "bl 

2m 



(6) 



with s = ± corresponding to |, |. We have neglected the 
real part of the Green's functions G R,A in the collision 
integral of the kinetic equation @ , since it only results 
in a small correction of order (rEp)^ 1 to the Rashba 
precession term. 

After solving the kinetic equation ©, we can find the 
spin density s(r, t) in terms of the Keldysh Green's func- 
tion g(r, t; p, e) (in the Wigner representation)^ 

<r,t) = ^ [ ^Tr[«T.g(r,t;p, e )] . (7) 



(2tt) 3 



where f3 = l/fc^T is the inverse temperature, and we 
will be interested in the T — > limit in the following. 
The kinetic equation can be simplified after defining the 
distribution function /(e, 9) at given position and time: 



g(e,p) = i{i(e,p),/M)} , 



(9) 



where we assume that at given energy e and momentum 
direction 9, the Keldysh Green's function g has the mo- 
mentum profile determined by the spectral function, i.e., 
it can be expressed as C|(e, 9)A^(p, e) + Cx(e, 9)Ai(p 7 e) 
in terms of some Hermitian matrix coefficients Cf,j.. 
Note that in equilibrium, f(e, 9) = — tanh(/3e/2)/2. 
Ansatz © allows us to integrate Eq. @ over the momen- 
tum (J pdp/2n) at fixed direction 9 and energy e. Then 
we can integrate over the energy e and look for a steady- 
state solution. To this end, we have to solve the following 
transport equation for the distribution function: 
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where f(9) — J def(e,9), ()g denotes angular averaging 
and j(y) = a(y)/vp- §e is the source term proportional 
to the electric field, 

§e = eE [— cos9 + 7(sin2# a x — cos 29 a y )] , (11) 

and we have to supply Eq. (jTOJ) with the "initial condi- 
tion," 



V = 0- f(8) = f + eElco S d. 



(12) 



Here, /o is an isotropic spin and charge imbalance (as- 
sumed to be linear in 7 as it vanishes for 7 = 0), which 
can in general be induced by the electric field in the bulk 
and leak towards the normal boundary, jo is assumed to 
be angle 9 independent, since any nonisotropic compo- 
nent of fo decays on the scale of the mean free path at 
y < 0, if we take «i = 0. The second term in Eq. (TT2"|) 
is the usual drift along the x axis. We can now integrate 
Eq. (10]) using Eqs. (JI]) and $12j). For the BC across the 
transition region < y < d which complements the diffu- 
sion equation (that is second order in spatial derivatives), 



we need to derive the expressions for the spin density s 
and d y s at y = d, i.e., after 7 is fully turned on, for a 
given spin density at y = (which can then be deter- 
mined self-consistently after solving the diffusion equa- 
tion). Working to the lowest nontrivial order in 7, we 
have to calculate s to the first order and d y s to the sec- 
ond order, as will be explained below. 



III. BOUNDARY SPIN ACCUMULATION 



To linear order in 7, Eq. (jTTJJ) simplifies tremendously 



to 



§' E +sme d y f+eEl^l sin20 fj e = j ((f) - /) , (13) 



where 



s' E = eE 



cos6> + ^ (sin 26» a x - cos26> a y ) . (14) 
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Equation (TT3)) was obtained by inserting the normal so- 
lution /jv = eEl cos 9 into terms in Eq. (|10p which have 
prefactors linear in 7. Averaging Eq. (|13|) and also 
Eq. (fl3l) multiplied by sin over 0, we get 

(f(y)) e =fo + eEl^-a y + (cos 26 f{y)) g . (15) 

We can evaluate the last term in this equation by notic- 
ing that if we fix and d y j in Eq. (fT5)l . then impu- 
rity scattering equilibrates (on the scale of the mean free 
path) the distribution to 

f{y)^(f{y)) g -l(s' E + eEl^s\n29^ , (16) 

where the last term should actually also be neglected in 
the spirit of the approximation (which disregards correc- 
tions of order l/d <C 1 with respect to the leading term 
linear in 7). Combining this with Eq. (|15[) . we get 

/(») = /o + eEl cos [1 - 7 (y)fo] . (17) 

This gives for the spin density (for the case ct\ = 0) 

s(y) = s(0) + ueETa(y)y , (18) 

where v — m/2n is the density of states per spin. This 
is just the position-dependent bulk result^ offset by s(0). 
The charge density n(y) is uniform. It may appear quite 
surprising (although we will provide a physical explana- 
tion in the next section) that the spin density follows the 
Rashba parameter profile which changes rapidly on the 
scale of the spin-diffusion length l s . 

Inserting the distribution function (|17p into the terms 
which already have prefactors linear in 7 in Eq. (|10p gives 
the equation for f(y,9) to the second order in the SO 
interaction, which can be solved giving at y = d: 

d y ii = (19) 

for the charge density and 

d y s = 2mapx x s — ImveETOt^L (20) 

for the spin density. Eq. ([2H)l is equivalent to the con- 
dition of vanishing normal (i.e., along the y axis) spin- 
current density at y = d (to the second order in a)&Z- 

j a -y = 0, (21) 

where a = x,y,z label spin-density components in the 
Cartesian coordinates. Combining Eq. (|20|) with the dif- 
fusion equation^ 

d t s a + {s a - SayVeEra) jr a = D\/ 2 s a 

+ 4i?FTaV • [z x a x z s z — (a • z) z x s x z] (22) 

for the region with constant a (here, r" 1 = Ty 1 = 
t z 1 /2 = 2A 2 r) allows to solve for the spin density at 



y > d induced by the small electric field E. We can, 
however, immediately see that since the in-plane bulk 
spin density s y results in the cancellation of the two terms 
contributing to d y s in Eq. (|2H)l . the uniform bulk solution 
all the way down to y — d trivially satisfies both the bulk 
diffusion equation $2$ and the BC ([22). Using Eq. JTSJ), 
we can extend the bulk solution at y = d down to y = 0, 
finding that s(0) = 0. The final result is thus that the 
nonequilibrium spin density induced by the electric field 
vanishes in the normal region y < 0, and for y > fol- 
lows the bulk value corresponding to the local Rashba 
parameter ct(y), as the latter is slowly turned on. 

We point out that if the term proportional to the elec- 
tric field in Eq. ([2U)l is disregarded 10 (which could physi- 
cally be justified, e.g., when the nonequilibrium spin den- 
sity is induced optically at the edge in the absence of 
electric field), then these BC would in general describe 
interconversion of the y and z components of the spin 
density, which decays into the bulk on the scale of l s ac- 
cording to the diffusion equation (|2"2"|) . This, however, 
gives a wrong result for the spin density generated by an 
electric field, as was previously noted in the case of a re- 
flecting Hall-bar edge by Bleibaum,— who obtained the 
same BC as Eq. (|20|) . in that case, using a very different 
method, and also previously postulated by Mal'shukov et 
al£ Finally, we note that from the structure of Eq. (|20[) . 
it should be clear why we calculated densities to linear 
order in a while gradients to quadratic order: The latter 
are governed by the spin precession in the Rashba field 
which is linear in a and the spin-generation term due to 
the electric field which is quadratic in a. The second- 
order in a BC for the normal components of the density 
gradients thus correspond to the first-order in a solution 
of the spin-diffusion equation. 

IV. BOUNDARY HALL EFFECT 

In this section, we offer a simple physical explana- 
tion for the spin-density jump across the interface, which 
will also allow us to extend our findings to more general 
boundary configurations. Because there is no net drift 
velocity in the y direction, the term proportional to a' in 
Eq. (TTJ) has much less effect than the term proportional 
to a. The primary role of a' is to cause relaxation of 
s y , which takes place only on the very long time scale 
t s or the corresponding length scale l s . In order to un- 
derstand variations on the much shorter length scale d, 
we may therefore set a' = in -ffso- Then, s y is con- 
served and a(r) is equivalent to a vector potential in the 
x direction, 

A x (r) = ±(m/e)a(r), (23) 

for electrons with a y — ±1, in addition to a potential 
independent of a y , which is cx a 2 . Thus, electrons with 
(Ty = ±1 feel an effective orbital magnetic field in the z 
direction, given by 

B z = z • (V x A) = T{m/e)d y a . (24) 
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Due to the Hall effect in B z , an electric field in the x 
direction leads to a chemical potential drop in the y di- 
rection, which has opposite sign for the two spin states, 
leading to a difference of populations, and hence a jump 
in s y , across the interval < y < d. The chemical po- 
tential jump may be calculated from the standard Drudc 
formula for the Hall resistance. Multiplying this by the 
density of states v, we find from Eq. 

s v (y) ~ e v (0) = veE f dy'T(y')d y a{y') , (25) 
Jo 

where r(y') is the transport scattering time at position y' . 
In the case where r is independent of y, the right-hand 
side of Eq. (|25p for y > d is just vcEt(olq — ai), which 
is precisely the difference in the bulk polarizations in the 
two media, far from the boundary between them. The so- 
lution of the coupled spin-diffusion equation (|22p is then 
simple: The polarization s y has one constant value in the 
region y < and another in the region y > d, with an 
effective discontinuity at y ~ 0. There are no gradients 
of s y on the scale of the spin-diffusion length on either 
side of the boundary, and hence the diffusion equation 
leads to no polarization out of the plane. This is true, in 
principle, even if there is a difference in the electron den- 
sity on the two sides of the interface, provided that the 
electron mobility, and hence r, is constant throughout. 

In general, however, density differences caused by 
charging a gate over one half of the system will lead to 
variations in r and hence a gradient of r in the transi- 
tion region < y < d. This will cause the discontinuity 
in s y given by Eq. (|25[) to deviate from the difference in 
bulk polarizations, forcing the existence of gradients in 
s y and nonzero values of s z , in a region of the size of the 
spin-diffusion length, on either side of y = 0. Of partic- 
ular interest is the case where a » on one side of the 
boundary, and the spin-diffusion length is especially long 
in that region. Then, if there is a gradient of electron 
mobility in the transition region, there can be a nonzero 
value of s z which extends far into the "normal region" 
where a w 0. The value, however, will depend on the 
position dependence of r as well as the behavior of a in 
the transition region. 



of an effective magnetic field B z very close to the bound- 
ary would have no physical effect: The generated Hall 
voltage is added to a potential which is in any case infi- 
nite on the insulating side of the interface. The solution 
of the diffusion equation ([22)1 for constant a at y > d 
with BC (|20p gives a spin polarization s y which is con- 
stant and the same as in the bulk, for points further than 
d from the boundary. We find that s z = in this case, 
which is true even when d — > 0, in agreement with Refs. 9 
and I 111. 

Finally, we note that in the case where the electric field 
is normal to a transmitting boundary, there should be no 
discontinuity in the spin polarization across the interface. 
If a had depended on x rather than y, in the discussion 
above, the vector potential A x would have been purely 
longitudinal, giving no physical effect. 

In summary, we have developed a method, starting 
from the quantum kinetic equation, for studying spin Hall 
effects in a system with a position-dependent Rashba cou- 
pling constant a, in the "dirty limit" At <C 1. We con- 
sidered systems where a depends only on y and changes 
discontinuously on the scale of the spin-diffusion length 
l s , from one constant value to another, with a uniform 
applied electric field E in the x direction, and we derived 
boundary conditions for the coupled spin and charge dif- 
fusion equations, which apply away from the discontinu- 
ity. In general, we find a discontinuity in the y compo- 
nent of the spin density s(r) at the boundary. In the 
case where the transport scattering time r is indepen- 
dent of position, this leads to values of s away from the 
boundary that are the same as would occur for a uniform 
system with the local value of a, namely, s v = veEra 
and s x = s z = 0. Thus, in contrast to Ref. 1 10. we find no 
lateral spin injection into a normal conductor and no out- 
of-plane spin density, when the electron mobility is uni- 
form. The related spin Hall current extraction proposal 
of Ref. |8| can therefore only work under the assumption 
of an inhomogeneous mobility close to the interface. We 
also find s z — at a reflecting boundary, in contrast 
to Ref. 10, but in agreement with Refs. and [Til On 
the other hand, for a transmitting boundary, if t varies 
across the transition region, one can get spin injection 
and s z 7^ 0. 



V. CONCLUSIONS AND DISCUSSIONS 

We have assumed so far that the scale d for variation 
of a is larger than the electron mean free path. In the 
case where r is a constant, however, our results should 
apply even to steplike variations in a, as long as we are 
not interested in the oscillating structure on the scale of 
the Fermi wavelength XpJ^ 

In the case of a perfectly reflecting boundary at y = 0, 
with current flow parallel to the boundary, the existence 
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